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Abstract 

I review recent progress in understanding the arrest and flow behaviour 
of colloidal glasses, based on mode coupling theory (MCT) and related ap- 
proaches. MCT has had notable recent successes in predicting the re-entrant 
arrest behaviour of colloids with short range attractions. Developments based 
upon it offer important steps towards calculating, from rational foundations 
in statistical mechanics, nonlinear flow parameters such as the yield stress 
of a colloidal glass. An important open question is why MCT works so well. 



1 Introduction 

1.1 Soft Matter 

This paper addresses issues of arrest and flow in soft condensed matter; sec |l| 
for a useful ensemble of background reading and Q for a good, experimentally 
motivated overview. We consider a system of N particles in the size range of 
nanometres (e.g. globular proteins) to microns (traditional colloids), suspended 
in a solvent, with total volume V. To a good enough approximation, at least for 
equilibrium properties, the solvent degrees of freedom can be integrated out to 
give an effective pairwise Hamiltonian H — J2i>j u ( r v) ( m an obvious notation). 
Equilibrium statistical mechanics, when applicable, is governed by the partition 
function Z — J e~xp(—f3H)T>[ri}; quantum effects play no role ||. 

In many colloidal materials the effective interaction u(r) comprises a hard 
core repulsion, operative at separation 2a with a the particle radius, combined 
with an attraction at larger distance. For simplicity one can imagine a square 
well potential of depth e and range £et, with £ < 1. Unlike atomic systems, to 
which colloidal ones are otherwise quite analagous, these parameters can be varied 
easily in experiment, essentially by varying the solvent conditions j^]. For example, 
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adding polymers to a colloidal system will mediate an entropic attraction between 
spheres whose range is comparable to the size of the polymer coils and whose depth 
is controlled by their concentration. For globular proteins the same tricks can be 
played with salt concentration and pH. 

The resulting equilibrium phase diagrams are well known, and depend on the 
range parameter £ and the attraction energy or temperature through the parameter 
f3e. For £ >0.2, at small j3e there is a phase separation from a colloidal fluid to a 
colloidal crystal. At higher /3e, a liquid phase intervenes; the fluid undergoes a gas- 
liquid separation at intermediate densities although the crystal is stable at higher 
ones. So far, this is just like the phase diagram of argon or a similarly classical 
atomic substance. (But of course solvent fills the space between the colloids, so 
the gas is not a real gas.) However, for smaller £ the liquid phase is missing: one 
has only one transition, from fluid to crystal. In fact, though, the liquid phase is 
still lurking beneath: it is metastablc. 

1.2 Arrest in Colloidal Fluids 

Colloidal fluids can be studied relatively easily by light scattering M . This allows 
one to measure the dynamic structure factor S(q, t) = (p(q, f)p(— q, t' + t))/N and 
also the static one, S(q) = S(q,0). Here p(r, t)) = X)i^( r «W — r ) — N/V; this is 
the real space particle density (with the mean value subtracted), and p(q, t) is its 
Fourier transform. For particles with hard-core repulsions, S(q) exhibits a peak 
at a value q* with q*a = 0(1). The dynamic structure factor S(q,t), at any q, 
decays monotonically from S(q) as t increases. In an ergodic colloidal fluid, S(q, t) 
decays to zero eventually: all particles can move, and the density fluctuations have 
a finite correlation time. In an arrested state, which is nonergodic, this is not true. 
Instead the limit S(q, oo) /S(q) = f(q) defines the nonergodicity parameter. (Note 
that this corresponds to the Edwards- Anderson order parameter in spin glasses.) 
The presence of nonzero f(q) signifies frozen-in density fluctuations. Although f(q) 
is strongly wavevector dependent, it is common to quote only f(q*) ||. 

Colloidal fluids are found to undergo nonergodicity transitions into two dif- 
ferent broad classes of arrested noncquilibrium state. One is the colloidal glass, in 
which arrest is caused by the imprisonment of each particle in a cage of neigh- 
bours. This occurs even for e = (i.e. hard spheres) at volume fractions above 
about <\> = 4ira 3 N/3V ~ 0.58. Such a system would, in equilibrium, be a crystal; 
but equilibrium can be delayed indefinitely once a glass forms (particularly if there 
is a slight spread in particle size a, which helps suppress nucleation). The noner- 
godicity parameter for the colloidal glass obeys f(q*) — 0.7. The second arrested 
state is called the colloidal gel. Unlike the repulsive glass, the arrest here is driven 
by attractive interactions, resulting in a bonded network structure. Such gels can 
be unambiguously found, for short range attractions, whenever (3e > 5 — 10. Hence 
it is not necessary that the local bonds are individually irreversible (this happens, 
effectively, at /3e Si 15 — 20); and when they are not, the arrest is collective, not local. 
It is found experimentally that for colloidal gels, /(g*)<;0.9, which is distinctly 
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different from the colloidal glass. The arrest line for gel formation slices across 
the equilibrium phase diagram (e.g., plotted on the ((f), j3e) plane), and, depending 
on £, parts of it lie within two phase regions. This, alongside any metastable gas- 
liquid phase boundary that is present, can lead to a lot of interesting kinetics ^, ^, 
in which various combinations of phase separation and gelation lead to complex 
microstructures and time evolutions. 



1.3 A Brief Primer on Mode Coupling Theory (MCT) 

This is not the place to explain MCT in detail. The most powerful form of the 
theory which is favoured by most of the true experts, remains somewhat ob- 
scure to many others. However, in a stripped down version (see e.g. ||, [l0|) the 
theory can be viewed as a fairly standard one-loop selfconsistent approach to an 
appropriate dynamical field theory. 

We take /3 = 1 and start from the Langevin equations = Fj + fj for 
independent particles of unit diffusivity (Dq = 1) subjected to external forces Fj. 
The noise force then obeys (fjfj) = l5ij. By standard manipulations one proceeds 
to a Smoluchowski equation $ = SI'S for the iV-particle distribution function 'J, 
with evolution operator is = Vj.(Vj — Fj). Now take the forces Fj to derive 
(via Fj = — Vji?) from an interaction Hamiltonian 



H =-^J <2 3 rd 3 r'p(r)p(r')c(|r - r'| 



(1) 



where Nc(q) = V[l — S(q) ]. This is a harmonic expansion in density fluctuations; 
c(q) is called the direct correlation function, and its form is fixed by requiring that 
S(q) be recovered in equilibrium. Neglected are solvent mediated dynamic forces 
(hydrodynamic couplings); these mean that in principle the Langevin equations 
for the particles should have correlated noise. Also neglected are anharmonic terms 
in H; to regain the correct higher order density correlators (beyond the two point 
correlator S(q)) in equilibrium, these terms would have to be put back. 
These assumptions give a Langevin equation for the density p(r): 

p = V 2 p + V(pV6H/Sp) + V.h (2) 

where h is a suitable noise (actually with a nontrivial density dependence pi). 
This equation is nonlinear, even with the harmonic choice of H. However, from it 
one can derive a hierarchy of equations of motion for correlators such as S(q,t), 
more conveniently expressed via $(q, t) = S(q,t)/S(q). Factoring arbitrarily the 
four-point correlators that arise in this hierarchy into products of two <&'s, one 
obtains a closed equation of motion for the two point correlator as 



Hq,t)+T(q) 



$(q,t)+ / m(q,t-t')$(q,t) 



(3) 
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where T(q) = q 2 /S(q) is an initial decay rate, and the memory function obeys 

m(q,t)=2V q>k *(k,i)$(k-q,t) (4) 

k 

with the vertex 

V «* = ^^S(q)S(k)S(\k q|)[q.kc(fc) + q.(k - q)c(|k - q|)] 2 (5) 

1.4 Dynamic Bifurcation 

The MCT equations exhibit a bifurcation that corresponds to a sudden arrest 
transition, upon smooth variation of either the density <f> or other parameters 
controlling the kernel in the harmonic hamiltonian H . This is best seen in the 
nonergodicity parameters f(q), which suddenly jump (for all q at once) from zero 
to nonzero values. Near this (on the ergodic side, which is always the direction 
MCT approaches from), $(q,t) develops interesting behaviour. Viewed as a func- 
tion of time, it decays onto a plateau of height f(q), stays there for a long time, and 
then finally decays again at very late times. The two decays are called (3 and a re- 
spectively. Upon crossing the bifurcation, the a relaxation time diverges smoothly 
with the parameters; upon crossing the locus of this divergence, f(q) = S(q, oo) 
therefore jumps discontinuously from zero to a value that is finite for all q. 

2 Critiques and Defence of MCT 

This mathematical structure means that S(q) fixes all of the dynamics (up to a 
scale factor which was the bare diffusion constant of one particle, here set to unity). 
As a result, the theory has significant predictive power: MCT gives, in terms of 
S(q) , the critical parameter values where arrest occurs; the power law exponents 
governing the a and j3 decays; the divergence of the a relaxation time; and f(q). 
This makes the theory falsifiable. Indeed, it has successfully been falsified. For 
example, there is no doubt that MCT gives the wrong density (f> g for the glass 
transition in hard spheres. On the other hand, the predictions (including that for 
4> g ) do all agree with experiments at about the 10 percent level, and often better. 

Attitudes among theorists to MCT for colloids therefore vary considerably. 
Some argue that the approximations made are uncontrolled (true) and that, even 
though there are no explicit adjustable parameters in the theory, its approxima- 
tions have been implicitly tuned to suit the problem at hand (normally taken 
to be the glass transition in hard spheres). These might be fair accusations in 
part, but if so they can also be levelled at much of equilibrium liquid state theory, 
where, for decades, ad-hoc closures (Percus-Yevick, hypernetted chain, etc.) have 
competed for survival by an essentially Darwinian process. Other theorists point 
to the already-falsified status of MCT, drawing attention to the misplaced </> s or, 
more interestingly, to other physical situations where exactly the same kind of 
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approximation leads to totally wrong predictions (e.g. spurious arrest in systems 
that are known to evolve smoothly towards a unique equilibrium state). 

Against all this must be weighed the continuous stream of experimental re- 
sults, many of them subtle but others not, that confirm the MCT predictions in 
surprising detail — especially if one makes one or two ad-hoc 'corrections' to adjust 
parameters like 4> g ]T^ ] . (There are also many simulation results which increasingly 
confirm the same picture @].) 

A striking recent success of MCT concerns systems with attractive interac- 
tions as well as hard-core repulsions, in a region of parameter space where not 
one but two arrest transitions are nearby (these are the transitions to a colloidal 
glass and to a colloidal gel). First, MCT unambiguously predicts that adding 
a weak, short range attraction to the hard sphere system should melt the glass. 
The mechanism appears to be a weakening of the cage of particles in which a given 
particle is trapped. (This is caused by members of that cage moving closer together 
under the attractive forces; gaps in the cage, allowing motion, are then more likely 
to appear.) Second, MCT predicts that adding more of the same attraction should 
mediate a second arrest, this time into a gel (f(q*) ~ 0.95). Third, MCT predicts 
that as parameters are varied, a higher order bifurcation point should be seen, 
when the re-entrant arrest line (implied by the above picture) crosses from being a 
smooth curve to being a cuspy one Jl3| . Although not every detail of this scenario 
is yet confirmed, there is clear experimental evidence of the predicted re-entrant 
behaviour connecting the glass and the gel arrest lines Q (see Figure [l]) and 
clear evidence of the proximity of the higher order bifurcation, which shows up as 



a characteristic logarithmic decay for t) 15 . The latter is also seen clearly in 
recent simulations [|l6| ; see Figure ^. 

It is worth emphasizing that these results were predicted by MCT before 
the experiments (or simulations) were actually begun. The success of MCT at 
unifying the glass and gel arrest transitions in attractive colloids does much to 
dispel fears of implicit parameter tuning during the earlier evolution of the theory. 
It reconfirms the generic success of the MCT approach, at around the 10 percent 
level, in describing interacting colloids. 

So, perhaps the time has come to stop criticising MCT for colloids on the 
grounds that it 'cannot be right'. Indeed it cannot: the problem is really to under- 
stand how it can do so well. To paraphrase Churchill |l7]]: Mode coupling theory 
is the worst theory of colloidal glasses - apart from all the others that have been 
tried from time to time. 



3 Shear Thinning 

In Refs.]l8t, M. Fuchs and the author developed a theory, along MCT lines, of col- 
loidal suspensions under flow. The work was intended mainly to address the case 
of repulsion-driven glasses, and to study the effect of imposed shear flow either on 
a glass, or on a fluid phase very near the glass transition. In either case, simpli- 
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Figure 1: Experimental test of MCT for colloids with short range attractions. The 
dashed line separates ergodic samples (open symbols) from nonergodic samples, 
identified as gels (squares), glasses (circles) and intermediate (+). The control 
parameter on the vertical axis is the concentration of added polymer; horizontal 
is the volume fraction. The solid line is the MCT prediction for the nonergodicity 
transition (shifted to give the correct <f> g in the absence of polymer). Courtesy K. 
Pham; see Poon et. al. |Q for individual discussions of the samples marked A-G. 

fications might be expected because the bare diffusion time To = a 2 /Do is small 
compared to the 'renormalized' one r = a 2 /D, which in fact diverges (essentially 
as the a relaxation time) as the glass transition is approached. If the imposed 
shear rate (which we assume steady) is 7, then for 770 « 1 < 77", one can hope 
that the details of the local dynamics are inessential and that universal features 
related to glass formation, should dominate. Note, however, that by continuing 
to use a quadratic H (Eq.Q), we will assume that, even under shear, the system 
remains 'close to equilibrium' in the sense that the density fluctuations that build 
up remain small enough for a harmonic approximation to be useful. 

The basic route followed in Ref.[[l8| is quite similar to that already laid out 
above for standard MCT. However, we assume that an imposed shear flow is 
present; obviously this changes the equations of motion. A key simplification is to 
neglect velocity fluctuations so that the imposed shear flow is locally identical to 
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Figure 2: Simulation test of MCT for Brownian spheres with short range attrac- 
tions. The correlator Q(q, t) is shown for various q, in a system thought to lie close 
to the higher order singularity at the glass-gel corner. Dotted lines are fits to the 
logarithmic time dependence predicted by MCT in this region. (From |fl6|.) 



the macroscopic one; this cannot be completely correct, but allows progress to be 
made. 
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For related earlier work see Refs. 
We again take (3 = 1, Dq = 1, and start from the Langevin equations = 
u + Fi + fj for independent particles of unit diffusivity subjected to external forces 
Fj and, now, an imposed flow velocity u(rj). We take this to be a simple shear 
flow with u(r) = jyx. The Smoluchowski equation ^ = is unchanged but the 
evolution operator is now f2 = ^ Vj.(Vj — Fj — u(rj)). We again take the forces 
Fi to derive from Eq.^, and with the same assumptions as before gain a Langevin 
equation for the density p(r): 



p + u.Wp 



w 



V{pVSH/Sp) + V.h 



(6) 



So far, the adaption to the equations to deal with imposed shearing is fairly 
trivial. The next stages are not. We assume an initial equilibrium state with <3>(i = 
0) oc Z , and switch on shearing at t = 0+. We define an advected correlator 



*(q,*) = (p(<l ) 0)p(-q(t),t))/S(q)N 



where q(i) = q + q.Ki with K the velocity gradient tensor, Kij 
definition of the correlator subtracts out the trivial part of the advection, which is 



(7) 

r y5 ix 5j y . This 



7 



merely to transport density fluctuations from place to place. The nontrivial part 
comes from the effect of this transport on their time evolution; the main effect 
(see e.g. ^0|) is to kill off fluctuations by moving their wavenumbers away from 
q* where restoring forces are weakest (hence the peak there in S(q)). Hence the 
fluctuations feel a stronger restoring force coming from H, and decay away more 
strongly. This feeds back, through the nonlinear term, onto the other fluctuations, 
including ones transverse to the flow and its gradient (i.e., with q along z) for 
which the trivial advection is absent. 

There follow a series of MCT-like manipulations which differ from those of 
the standard approach because they explicitly deal with the switchon of the flow 
at t = 0+ and integrate through the transient response to obtain the steady state 
correlators, under shear, as t — ► oo. There is no integration through transients 
in standard MCT; on works directly with steady-state quantities. (In practice 
also, the following results were obtained in Refs. [|l8) using a projection operator 
formalism which differs in detail from the version of MCT outlined above.) Despite 
all this, the structure of the resulting equations is remarkably similar to Eqs. pH: 



*(q,t)+r(q,t) 



$(q,t)+ / m(q,M')$(q,0 
Jo 



(8) 



with T(q) replaced by a time dependent, anisotropic quantity: 

T(q, t)S{q) = q 2 + q x q y "ft + (q x q y jt + qlj 2 t 2 )S(q) - q x q y S{q)/q (9) 

The memory kernel is no longer a function of the time interval t — tf but depends 
on both arguments separately 



m 



; q , t, t') = Y, V(l, k > t, 0«(k, - q, t) (10) 



through a time-dependent vertex V, too long to write down here (Taj. 

Using a nonequilibrium Kubo-type relationship, one can also obtain an ex- 
pression for the steady state viscosity r\ = cr{ A f)/'j where (7(7) is the shear stress 
as a function of shear rate. The viscosity is expressed as an integral of the form 

/>OC 

n= / dty>(k,t)$ a (k,t) (11) 

Jo k 

where the function F may be found in Ref. O] . 



3.1 Results 

The above calculations give several interesting results. First, any nonzero shear 
rate, however small, restores ergodicity for all wavevectors (including ones which 
are transverse to the flow and do not undergo direct advection). This is important, 
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since it is the absence of ergodicity that normally prevents MCT-like theories being 
used inside the glass phase, at T < T g or <f> > <f) g . Here we may use the theory in 
that region, so long as the shear rate is finite. 

In the liquid phase (jfi < 4> g ) the resulting flow curve 17(7) shows shear thin- 
ning at 7T<;1, that is, when the shearing becomes significant on the timescale of 
the slow relaxations. This is basically as expected. Less obviously, throughout the 
glass, one finds that the limit a(j — » 0+) = cry is nonzero. This quantity is called 
the yield stress and represents the minimum stress that needs to be applied before 
the system will respond with a steady-state flow. (For lower stresses, various forms 
of creep are possible, but the flow rate vanishes in steady state.) 

The prediction of a yield stress in colloidal glasses is significant, because 
glasses, operationally speaking, are normally defined by the divergence of the vis- 
cosity. However, it is quite possible for the viscosity to diverge without there being 
a yield stress, for example in 'power law fluids' where ~ j p with < p < 1. 
Indeed, a recent model of 'soft glassy materials' (designed for foams, emulsions, 
etc., and based on the trap model of Bouchaud |gl[) gives a well-developed power 
law fluid region above T g in which the viscosity is infinite but the static shear 
modulus zero p^ |. This does not happen in the present calculation, where the 
yield stress jumps discontinuously from zero to a nonzero value, cr y , at 4> g . The 
existence of a yield stress seems to be in line with most experimental data on the 
flow of colloidal glasses, although one must warn that operational definitions of 
the yield stress do vary across the literature ^3|. Ours is defined as the limiting 
stress achieved in a sequence of experiments at ever decreasing 7, ensuring that a 
steady state is reached for each shear rate before moving onto the next one. The 
latter requirement may not be practically achievable since the equilibration time 
could diverge (certainly one would expect to have to wait at least for times t such 
that 7i>I). But unless the flow curve has unexpected structure at small shear 
rates, the required extrapolation can presumably be made. 

The existence of a yield stress everywhere within the glass phase follows from 
the structure of the MCT-inspired calculations outlined above and detailed in Ref. 
[ fl8| . However, to calculate an actual value for oy requires further approximations; 
these avenues are pursued in Ref. jl8|. Quantitative results for one such approx- 
imation, called the 'isotropically sheared hard sphere model' (ISHSM), are given 
in Figure [|. Such approximations can also give values for the flow curve exponent 
after the onset of yield (a — cry ~ 7 P with p ~ 0.15), and for the growth of the 
yield stress beyond the glass transition (cry — a Y ~ {4> — (fig) 1 / 2 ). 

3.2 Schematic MCT models 

It has long been known that the key mathematical structure behind (conventional, 
unshcared) MCT can be captured by low-dimensional schematic models in which 
the full q dependence is suppressed. In other words, one chooses a single mode, 
with a representative wavevector around the peak of the static structure factor, and 
writes mode coupling equations for this mode in isolation. At a phenomenological 
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Figure 3: Flow curves for the ISHSM. The shear stress is given in units of kT/d 3 
with d the particle diameter. The parameter Peo is jd 2 /Dq with Dq the bare 
diffusion constant. The curves marked with circles are within the fluid phase (with 
dashed asymptotes showing the Newtonian limit); the curves marked with squares 
are in the glass. The diamonds denote the critical case. Each curved is marked by 
its distance from the glass transition, <fi — (f)g. The yield stress at the glass transition 
<jy (here called af) is indicated on the left. 



level, one can capture the physics similarly even with shearing present (despite the 
more complicated vectorial structure that in reality this implies). Specifically one 
can define the F^ 2 model — the sheared extension of a well known static model, 
F 12 — via 



*(*) 



m{t - t')$>(t')dt' 



= 



with memory function (schematically incorporating shear) 
m(t) = [vi$(t) + w 2 $ 2 (i)]/(l + 7 2 i 2 ) 



(12) 



(13) 



The vertex parameters are smooth functions of the volume fraction (j) (and 
any interactions). To calculate flow curves, etc., one also needs a schematic form 
of Eq.0; here we take the first moment of the correlator to fix the time scale for 
stress relaxation (which is, in suitable units, simply the viscosity): 



(14) 
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(Note that a different choice, e.g. with $(i) 2 in this equation to closer resemble 
Eq.jll], would yield quite similar results.) This schematic model gives very similar 
results to the ISHSM, with a - a Y ~ 7 016 and a Y - a Y ~ O - 4> g ) 1/2 @. The 
qualitative reproducibility of these results within different types of approximation 
scheme is reassuring. 



4 Shear Thickening and Jamming 

The calculations described above predict, generically, shear thinning behaviour: 
advection kills fluctuations, reducing the a relaxation time, which causes the sys- 
tem to flow more easily at higher stresses. However, in some colloidal systems, the 
reverse occurs. This is shear thickening, and gives a flow curve cr(j) with upward 
curvature. In extreme cases, an essentially vertical portion of the curve is reported 
[ p4[ . One interpretation of the latter scenario (called 'discontinuous shear thick- 
ening') is that the underlying flow curve is actually S-shaped. Since any part of 
the curve with negative slope is mechanically unstable (a small increase in the 
local shear rate would cause an acceleration with positive feedback) , this allows a 
hysteresis cycle in which, at least according to the simplest models, discontinuous 
vertical jumps on the curve bypass the unstable section (see Figure ||). 

If this viewpoint is adopted, there seems to be nothing to prevent the upper, 
re-entrant part of the curve to extend right back to the vertical axis (see Figure 
|J) in which case there is zero steady-state flow within a certain interval of stress. 
The system has both an upper and a lower yield stress delimiting this region. (If 
it is nonergodic at rest, it could also have a regular yield stress on the lower part 
of the curve near the origin - we ignore this here.) This case has been called 'full 
jamming' p5| . Although mostly a theoretical speculation, one or two experimental 
reports of this kind of behaviour have appeared in the literature recently [ p6| . 

The above discussion suggests that shear thickening and full jamming might 
be viewed as a stress-induced glass transition of some sort. If so, it is natural to 
ask whether this idea can be accommodated within an MCT-like approach. Since 
the analysis of Ref. (ll| gives only shear thinning, this is far from obvious. In 
particular, a stress-induced glass transition would require the vertex V to 'see' the 
stress; this might require one to go beyond harmonic order in the density, that is, 
it might require improvement to Eq.|l|. Indeed, since it is thought that jamming 
arises by the growth of chainlike arrangements of strong local compressive contacts 
[ p7| , it is reasonable that correlators beyond second order in density should enter. 
In very recent work, an ad-hoc schematic model along the lines of Eqsjl2-14 



has been developed to address shear thickening. The vertex is ascribed explicit 
dependence not only on 7 (as suggested by the shear thinning calculations of Ref. 
[fl8|) but also on the shear stress a. It is found that a vertex which is monotoni- 
cally decreasing with 7 but monotonically increasing with a can indeed result in 
shear thickening, and, under some conditions, in full jamming p8[ . This work is 
preliminary, but interesting in that it suggests how new physics (beyond two-point 
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Figure 4: Three possible flow curves for a shear thickening material. The monotonic 
curve corresponds to continuous shear thickening. The remaining two curves are 
S-shaped; one expects, on increasing the shear rate, the stress to jump from the 
lower to upper branch at (or before) the vertical dashed line shown in each case. 
One curve shows the full jamming scenario: the existence of an interval of stress, 
here between 0.45 and 0.63, within which the flow rate is zero, even in a system 
ergodic at rest. (Stress and strain rate units are arbitrary.) 



correlations) may need to be added to MCT before the full range of observed col- 
loidal flow behaviour is properly described. Of course, even for systems at rest, it 
is known that some important physics is missing from MCT, in particular various 
kinds of 'activated dynamics' in which the system can move exponentially slowly 
despite being in a region of phase space where, according to MCT, it cannot move 
at all (see e.g. |Kj). Jamming seems very different from this: so perhaps there are 
more things missing from MCT than just activated processes. 

5 Conclusion 

Mode Coupling Theory (MCT) has had important recent successes, such as pre- 
dicting, in advance of experiment, the re-entrant glass/gel nonergodicity curves 
that arise in colloidal systems with short range attractions |u| [l4|] . 

Theoretical developments directly inspired by MCT now offer a promising 
framework for calculating the nonlinear flow behaviour of colloidal glasses and 
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glassy liquids fll8| . In fact, this is the only quantitative framework currently in 
prospect for the rational prediction of yield behaviour and nonlinear rhcology in 
this or any other class of nonergodic soft materials. (Other work on the rheology 
of glasses [^2| |2^] does not, as yet, offer quantitative prediction of experimental 
quantities.) While promising, many things are missing from our approach: velocity 
fluctuations, hydrodynamic forces, anharmonicity in H (possibly implicated in 
shear-thickening) etc., are all ignored. Further hard theory work is needed here. 

The eventual goal lies beyond the steady state flow curve (shear stress as 
a function of shear rate, (7(7)) discussed in this article. Ideally we would like a 
full constitutive equation that relates the stress tensor at a given time to the 
preceding deformation history (or vice versa). To do this, starting from statistical 
mechanics, is difficult for any class of material; its achievement for the case of 
entangled polymers j30), which are ergodic, was a highlight of theoretical physics 
in the late 20th century. To obtain a well-founded constitutive equation for a 
significant class of nonergodic soft materials is a worthy goal for the early 21st. 

As a framework for theoretical prediction, MCT is easy to criticize, but much 
harder to improve. The key question that should exercise the minds of theorists 
is no longer 'is it correct' (it clearly is not, in a technical sense) but 'why does it 
work so well' ? Only by engaging with this question are we are likely to improve 
our understanding of nonergodic colloidal materials and, in the longer run, come 
up with something better. 

Acknowledgements: I am deeply indebted to Matthias Fuchs, much of whose 
work is reviewed above, for introducing me to MCT. 
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